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ABSTRACT 

We study how the sizes and radial profiles of galaxies vary with wavelength, by fitting 
Sersic functions simultaneously to imaging in nine optical and near-infrared bands. To 
quantify the wavelength dependence of effective radius we use the ratio, TZ, of mea¬ 
surements in two restframe b ands. The dependence of Sersic index on wavelength, M, 
is computed correspondingly. Vulcani et al. (2014) have demonstrated that different 
galaxy populations present sharply contrasting behaviour in terms of TZ and M. Here 
we study the luminosity dependence of this result. We find that at higher luminosi¬ 
ties, early-type galaxies display a more substantial decrease in effective radius with 
wavelength, whereas late-types present a more pronounced increase in Sersic index. 
The structural contrast between types thus increases with luminosity. 

By considering samples at different redshifts, we demonstrate that lower data 
quality reduces the apparent difference between the main galaxy populations. However, 
our conclusions remain robust to this effect. 

We show that accounting for different redshift and luminosity selections partly 
reconciles the size variation measured by [Vulcani et al.| with the weaker trends found 
by other recent studies. Dividing galaxies by visual morphology confirms the behaviour 
inferred using morphological proxies, although the sample size is greatly reduced. 

Finally, we demonstrate that varying dust opacity and disc inclination can account 
for features of the joint distribution of TZ and N for late-type galaxies. However, dust 
does not appear to explain the highest values of TZ and Af. The bulge-disc nature of 
galaxies must also contribute to the wavelength-dependence of their structure. 

Key words: galaxies: general - galaxies: structure - galaxies: fundamental parame¬ 
ters - galaxies: formation 
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1 INTRODUCTION 

The assembly history of a galaxy is written in its stellar pop¬ 
ulations. Their colours (and spectral features) reflect the 
time and conditions of star formation, while their spatial 
(and kinematic) distributions provide insight into the struc¬ 
tures in which the stars formed and their subsequent redistri¬ 
bution by mergers and secular processes. These observables 
are all coupled; studying the detailed relationships between 
them helps us identify typical scenarios and unusual cases. 
We can therefore hope to understand more about the de¬ 
velopmental histories of galaxies and their internal compo¬ 
nents by studying the wavelength dependence of their spatial 
structure. 

Galaxies can be split visually into two main morpholog¬ 
ical types: elliptical and spiral. The primary distinguishing 
feature is the presence of a disc component. Lenticular galax¬ 
ies are commonly grouped with ellipticals, as early-types, al¬ 
though their disc-component suggests they are more closely 
related to spirals ( van den Bergh|1976 1. To complicate mat¬ 
ters, even ellipticals sometimes appear to have a somewhat 
‘disky’ appearance i Kormendy fc Bender|19^ l or kinemat¬ 
ics ( Cappellari et al.||2011 l. 

The radial luminosity profile of a galaxy, or its compo¬ 
nents, is often modelled by a Sersic profile. This function 
provides a good description of how the surface brightness 
of many galaxies varies with radius ( Sersic|1963 Graham & 
Driver|2005 


Elliptical galaxies are typically well-described by a sin¬ 
gle Sersic profile. The index is often taken to be fixed at 
n = 4 ( |de Vaucouleurs||1948p , although it has been shown 
to vary from n < 2 at low masses to n > 4 for high mass 
ellipticals (Hodge|1971 Gaon, Capaccioli fc D’Onofrio|19^ 


|Graham|26l3[). 

Spiral galaxies, on the other hand, generally require a 
combination of two different Sersic profiles. The discs of spi¬ 
rals (and lenticulars) can usually be accurately described by 
an exponential (n = 1) prohle (e.g., Kormend^|1977 Allen 
et al.|2006 1. Bulges are often described by a de Vaucouleurs 


(n = 4) profile. As with elliptical galaxies, this is usually not 
accurate ([Andredakis, Peletier fc Balcells||1995 Graham & 


Prieto|2001[|de Souza, Gadotti fc dos Anjos|2004 l , so recent 

studies (including this one) tend to adopt a general Sersic 
profile, where n is not fixed. Some studies also account for 
the presence of a bar, commonly with a third Sersic profile 
| |Gadotti|2011[ ). 

Ideally we would be able to decompose galaxies into 
all their constituent structures. However, with limited phys¬ 
ical resolution and signal-to-noise, even bulge-disc decom¬ 
positions can prove difficult to perform robustly. Although 
many galaxies comprise multiple components, fitting with a 
single Sersic profile can provide valuable insights into their 
structure, providing sizes, axial ratios, rough morphologi¬ 


cal classihcations and total photometry (e.g., Kelvin et al. 


2012 ). 


Fitting Sersic models to galaxy images in different wave¬ 
bands allows one to explore the dependence of galaxy struc¬ 
ture on wavelength for large samples in a homogeneous man¬ 
ner. Variation in a galaxy’s surface brightness profile param¬ 
eters with wavelength is qualitatively equivalent to a radial 
trend in the colour of the galaxy. Such ‘colour gradients’ 
have been well-studied, usually by measuring the surface 


brightness in a series of elliptical annuli. This approach is 
usually limited by the seeing, and hence confined to rela¬ 
tively extended objects. The result is a colour profile for the 
galaxy, usually summarised as the overall change in colour 
over some (logarithmic) radial range (e.g., |Goudfrooij et al.| 


1994 deJong|1996| Saglia et al.|2000 den Brok et al.|2011 1. 

The wavelength-dependence of Sersic model parameters 
provides an alternative description. This is less detailed than 
a full colour profile, though more robustly determined in 
typical survey imaging. In addition, as the PSF can be easily 
accounted for when fitting the model, the results are not 
biased by seeing. 

Gonsidering the wavelength-dependence of Sersic pa¬ 
rameters is more efficient than dealing with empirical colour 
profiles, while providing significantly more general informa¬ 
tion than quoting a colour gradient over a single (or several) 
radial ranges. Furthermore, the ubiquity of the Sersic profile 
in describing galaxy profiles, and the correlations between 
its parameters and other galaxy properties, suggests that 
the profile and its parameters have a physical significance. 
Expressing colour gradients in terms of the wavelength- 
dependence of effective (half-light) radius and Sersic index 
naturally separates the proportions associated with changes 
in size and profile shape, which may have different physical 
drivers. However, the price of all this is the assumption of 
a particular functional form for the surface brightness pro¬ 
file. When this assumption is inappropriate (e.g., in merging 
galaxies, multi-component systems or potentially very high 
redshift), it may result in biases or misinterpretations. 

Some studies have chosen a compromise between per¬ 
forming a parametric fit and directly measuring a colour 


profile. La Barbera & de Garvalho (20091 fit independent 


Sersic models in various wavebands and then use elliptical 
annuli to measure colour profiles and gradients on the model 
images. This removes the effect of the PSF and reduces noise 
in the measurements. 

Several groups have examined the trends of Sersic model 
parameters with wavelength. For a sample of bright ellipti¬ 
cal galaxies, EE Barbera et al.| ( |2010b[ ) find that the mean 
effective radius decreases significantly with increasing wave¬ 
length from g to H, but that there is little variation in Sersic 
index. [Kelvin et al.| (|2012[) find similar behaviour in data 


from the Galaxy And Mass Assembly survey (GAMA; Liske 
et al.|20r5 I. They also consider late-types, finding a substan¬ 


tial increase in Sersic index, and decrease in effective radius, 
across the same wavelength range. Metallicity gradients and 
dust attenuation are proposed as reasons for these trends. 


Lange et al. (2015 I have recently extended the work of 


Kelvin et al. (2012 I, replacing the GAMA infrared UKIDSS 
LAS imaging ( [Lawrence et aL]|2007P with deeper VISTA¬ 
VIKING data (Edge et al. 2013). They also separate the 


galaxy population into early- and late-types using a variety 
of dividers, including visual classification. [Lange et ah] find 
that their results are not sensitive to chosen divider, but 
show that the bimodality of galaxy structure becomes less 
distinct at lower stellar masses. They confirm that galax¬ 
ies appear more compact at redder wavelengths, although 
this appears less dramatic than in the previously mentioned 
studies. [Lange et al.j argue that these structural variations 
with wavelength may arise from the two-component nature 
of many galaxies, in which the bulge is observed in the red¬ 
der wavebands whilst the disc is observed in the bluer wave- 
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bands (also see Driver et al.pOOTa l. However, they do not 
examine the dependence of Sersic index on waveband. 

The above studies fit the image in each filter-band in¬ 
dependently. This model freedom can be seen as an advan¬ 
tage. However, such an approach does not ntilise the ex¬ 
pected strong correlations between structural parameters 
at neighbouring wavelengths. In data with low signal-to- 
noise or poor resolution, using these physical expectations 
can improve the reliability and precision of structural mea¬ 
surements. Furthermore, even in high-quality data, utilising 
available colour information can improve the performance 
of decompositions. To address these issues, the MegaMorph 
project ( [Bamford et al.||2012p has developed a technique 
which fits a single, wavelength-dependent model simultane¬ 
ously to a set of imaging in different filter-bands. These de¬ 
velopments have been implemented in an extended version of 
GALFIT (Peng et al. 2002| 20101 named galfitm (Haufiler 
et al. 2013| hereafter H13; [Vika et al. 2013 20141. GALA- 
PAGOS ^Barden et al.||2012| , a software package dedicated 
to running GALFIT on large surveys, has been similarly ex¬ 
tended to become galapagos- 2, which uses galfitm (H13). 

Using MegaMorph multi-band techniques to perform 


single-Sersic fits, Vulcani et al. (2014, hereafter V14) stud¬ 
ied the wavelength dependence of galaxy effective radius 
and Sersic index. They confirmed the trends described 
above, and demonstrated that galaxies with contrasting to¬ 
tal colours and Sersic indices display strikingly different be¬ 
haviour in terms of the wavelength dependence of their ef¬ 
fective radii and Sersic indices. 


V14 speculate that the behaviour seen for early-types 
may result from the superposition of red compact structures 
and more-extended, bluer structures, all of which possess 
the n ~ 4 profile of systems supported by random motions. 
This is similar to the picture presented by |Huang et al.| 
(20131. The behaviour of late-types is attributed to their 
two-component nature, as mentioned above. 


In this paper we aim to gain a greater confidence in, 
and deeper understanding of, the trends of galaxy struc¬ 
tural parameters with wavelength as measured by the var¬ 
ious studies described above. This work will advise our ap¬ 
plication of multi-band structural fitting to the full set of 
high-quality imaging currently being assembled by GAMA, 
and our subsequent analysis. We build on the results of V14, 
first ensuring they are robust to redshift effects. 

For the first time, we examine the wavelength- 
dependence of galaxy structure as a function of luminosity. 
We also show that the measured trends persist for selec¬ 
tions based on visually classifed morphology, rather than 
morphological proxies. Next, we investigate the role of dust 
in driving the wavelength dependence of late-type galaxy 
structure. We conclude by discussing the various potential 
factors responsible for the variation of galaxy structural pa¬ 
rameters with wavelength. By taking the novel approach of 
considering the joint wavelength dependence of Sersic index 
and effective radius, we can begin to decouple the relation¬ 
ships between concentration, size and colour. 


The analysis has been carried out using a cosmology 
with {Q,m,Q,K,h) = (0.3,0.7,0.7) and AB magnitudes. 


2 DATA 


2.1 Parent sample and structural parameters 

The galaxy measurements used in this paper have previously 
been presented in H13, and a subsample of these were stud¬ 
ied in V14. A detailed description of the selection criteria, 
robustness of fits and properties of the sample can be found 
in those papers; here, we only give a brief overview. 

Our data is taken from the extended Galaxy And Mass 
Assembly survey (GAMA H) (Liske et al, submitted), which 
is the largest homogeneous, multi-wavelength dataset cur¬ 
rently available for the low-redshift Universe. The imaging 


data GAMA has assembled, from both SDSS (York et al. 
|2000p and UKIDSS ( [Lawrence et al.|2007| ), provides a consis¬ 
tent set of pixel-registered multi-wavelength data, covering 
the (ugriz) optical bands and the (YJHK) near-infrared 
(NIR) bands, respectively ( [Hill et al.||201l| . It has been 
demonstrated that these data have sufficient depth and res¬ 
olution to allow Sersic profiles to be fit to large samples 
( [Kelvin et al.|20T^ H13). 

As in V14, the data we use in this paper are limited 
to the G09 region of GAMA. We therefore only utilise ap¬ 
proximately one-fifth of the area covered by GAMA H. For 
efficiency, we chose to focus on a single region for our initial 
exploration of multi-band fitting techniques with GAMA. 
The application of MegaMorph methods to the full GAMA 
dataset, with UKIDSS-LAS replaced by VISTA-VIKING 
data, are ongoing. 

We use GALAPAGOS-2 with galfitm to fit a single 
wavelength-dependent 2D model to all images of a galaxy 
simultaneously, galfitm allows each parameter in the stan¬ 
dard Sersic light profile to vary as a function of wavelength; 
the user may choose the degree of smoothness they re¬ 
quire. In this work, galaxy magnitudes are allowed to vary 
freely, while Sersic index and effective radius are modelled as 
quadratic functions of wavelength. All other parameters are 
not permitted to vary with wavelength, galfitm accounts 
for seeing by convolving the model with a PSF before com¬ 
paring to the data. We utilise the PSFs produced by [Kelvin[ 

[er^ ( (M^ . 

Restframe values of Sersic index and effective radius 
for each wavelength have been obtained from the polynomi¬ 
als returned by galfitm. For the magnitudes, k-corrections 
have been performed using InterRest ([Taylor et al.[[200^. 


2.2 Sample selection 

We take the GALAPAGOS-2 output catalogue, however not 
every fit was necessarily accurate or meaningful, so a number 
of cleaning criteria were applied to eliminate galaxies for 
which the fit has violated our criteria in one or more bands. 
These criteria are identical to those used in H13 and V14, 
with the exception of Sersic index which is limited to 0.201 < 
n < 7.75. We have reduced the upper limit on Sersic index 
to avoid including suspected poor fits, with final parameters 
close to the fitting constraint (objects with very high n are 
typically poor fits, often to point sources, and are therefore 
removed). 

In this paper we use two different volume-limited sam¬ 
ples to avoid selection effects, which are illustrated in Fig.[^ 
For studying variation in galaxy properties with redshift, a 
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z 

Figure 1. Absolute r-band magnitude versus redshift for our 
parent sample, with our volume-limited selection boxes overlaid. 
The curve line indicates the primary apparent magnitude limit of 
the GAMA II redshift survey, r < 19.8. This corresponds to an 
absolute magnitude of Mr = —21.2 at 2 = 0.3 and Mr = —19.48 
at 2 = 0.15. 




Ur 

< 2.5 

Tlr 

> 2.5 


Colour 

no. 

% 

no. 

% 

Blue: 

(u — r) < 1.6 

1380 

9.9 

322 

2.3 

Green: 

1.6 < (tx — r) < 2.1 

3311 

23.9 

1047 

7.5 

Red: 

2.1 < {u — r) 

2523 

18.2 

5288 

38.1 


Table 1. Number count and fraction of the z < 0.3, Mr < —21.2 
volume-limited sample (total = 13871 galaxies), for different com¬ 
binations of colour and Sersic index. 

first volume-limited sample of 13,871 galaxies is taken with 
2 < 0.3, Mr < —21.2, in line with V14. Within this red¬ 
shift limit the restframe iT-band can be interpolated from 
measurements in the observed u-K wavelength range. 

For studying variation in galaxy properties with abso¬ 
lute magnitude, Mr, a second volume-limited sample (z < 
0.15, Mr < —19.48) of 5331 galaxies is used, allowing galax¬ 
ies to be considered over a wider range of absolute magni¬ 
tude. 

Note that the cleaning described above removes only 
0.42 / 0.39 per cent of our V14 / low-redshift samples. 

One of the aims of this project is to study how galaxy 
properties vary as a function of wavelength. It is expected 
that different types of galaxies will behave in different ways, 
so we subdivide our sample (as in V14) by both colour and 
Sersic index. The number of galaxies in each subsample are 
given in Tables [^&[^ for both volume-limited samples, after 
basic cleaning has been applied. 

Initially the samples are divided into ‘red’ and ‘blue’ at 
{u — r) = 2.1. We also want to separate the bluest galaxies, 
which may contain starbursts, so we further divide this ‘blue’ 




Ur 

< 2.5 

rir > 2.5 


Colour 

no. 

% 

no. % 

Blue: 

(u — r) < 1.6 

1323 

24.8 

134 2.5 

Green: 

1.6 < (tx — r) < 2.1 

1417 

26.6 

245 4.6 

Red: 

2.1 < (u — r) 

994 

18.6 

1218 22.8 

Table 2. 

Number count and fraction 

of the 2 

< 0.15, Mr < 


— 19.48 volume-limited sample (total = 5331 galaxies), for differ¬ 
ent combinations of colour and Sersic index. 
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Figure 2. (u-r) rest frame colour vs. Ur for the galaxies in our two 
samples. Lines illustrate the cuts we apply to divide the galaxies 
by colour and Sersic index. 

sample at (u — r) = 1.6 into ‘green’ and ‘blue’. Note that our 
‘green’ sample corresponds to the main population of star¬ 
forming galaxies, not the green valley. We have confirmed 
that altering these colour cuts does not affect our results. 
The galaxy sample is also split by Sersic index, in an attempt 
to separate discy galaxies from ellipticals, at = 2.5. These 
divisions can be seen in Fig. 


3 RESULTS 

V14 examined the variation of n and with wavelength 
in order to reveal trends of internal structure, and therefore 
formation history, for galaxies in different subsamples. For 
consistency we continue to use the notation Mg = n{x)/n{g) 
and TZg = Rc{x)/Rc{g), which express the variation for a 
given waveband, x, versus that in the g-band. For simplicity, 
we will only consider the ratio of Jf-band to g-band in this 
paper, as V14 shows that there is a consistent trend with 
wavelength and this pair has the longest robust wavelength 
baseline. We omit the band labels from M and TZ when 
discussing their general behaviour. 

If we convert our M and TZ trends to colour gradients 
we hnd good agreement with previous work, e.g. red, high- 
n galaxies show identical colour gradients to the passive, 
‘early-type’ sample studied in |La Barbera et al.| ( |2010a| ). We 
also see the same behaviour of colour gradient with luminos¬ 
ity. However, here we focus on M and TZ, which allows us 
to decouple variations in concentration and size with wave- 
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length, which are combined in a measurement of colour gra¬ 
dient alone. 

V14 concluded that, regardless of overall u — r colour, 
low-n galaxies display a very different behaviour to high- 
n galaxies, implying that both A/” and TZ depend strongly 
on Sersic index. They suggested that low-n discy (two- 
component) galaxies show a large change in Sersic index 
with wavelength because in bluer wavebands the Sersic index 
of the disc is being measured, whereas in redder wavebands 
the Sersic index of the bulge dominates. High-n galaxies are, 
on the other hand, likely to be bulge-dominated (closer to 
one-component), so there is little change in Sersic index with 
wavelength. There is, however, a large change in Ra with 
wavelength for high-n galaxies, which indicates that ellipti¬ 
cal galaxies contain a number of different pressure-supported 
stellar populations with different extent, possibly resulting 
from multiple merging events throughout a galaxy’s lifetime. 

3.1 Redshift effects on galaxy structure 

The sample used in V14 extends to Mr < —21.2, « < 0.3, 
so much of the sample is faint and of small angular size. 
This selection was partly chosen to demonstrate the power 
of multi-band fitting in this challenging regime. V14 showed 
that the measured N and TZ correspond to differences in the 
visual appearance of the galaxies. However, before exploring 
the luminosity dependence of these trends, we will further 
test the resilience of our and TZf measurements versus 
redshift. 

We split our bright volume-limited sample [z < 0.3, 
Mr < —21.2) into three redshift bins: 0.0 ^ z < 0.15, 
0.15 ^ 2 < 0.25 and 0.25 ^ 2 : < 0.3. These bins were cho¬ 
sen to span the redshift range of our volume-limited sam¬ 
ple, while ensuring a similar number of galaxies in each bin 
to permit meaningful comparisons (see Table [^. In Fig. 
we show the wavelength dependence of effective radius, for 
galaxies split by u — r colour, Ur, and redshift. Galaxies of 
all colours and Sersic indices are typically found to have 
TZf < 1, indicating they are smaller at redder wavelengths, 
while TZf = 1 (i.e. log TZf = 0) corresponds to no variation 
in size with wavelength. Here, TZf does not appear to change 
substantially with redshift in any subsample. Kolmogorov- 
Smirnov (KS) tests do indicate some significant differences 
between the TZf distributions in the different redshift sam¬ 
ples. However, these differences are generally small, particu¬ 
larly in comparison with the width of each distribution and 
the offsets between the low- and high-n subsamples. To de¬ 
termine whether an offset between subsamples can be con¬ 
sidered ‘large’ or ‘small’ we sum the standard deviations of 
the widest and narrowest distributions in quadrature. We 
then find the difference in the median value of log{TZf) in 
the highest and lowest redshift bins, as a fraction of the 
summed standard deviation. Here, there average offset is 
15.5% of the distribution widths, which can be considered 
small. 

What trends are present, indicate that at lower redshift 
(higher S/N and better resolution) we measure TZf closer 
to unity for low-n galaxies. The departure of TZf from unity 
may therefore be slightly overestimated in V14. High-n sys¬ 
tems appear unaffected, with perhaps the exception of blue, 
high-n galaxies. There are, however, very few galaxies in this 
subsample as blue spheroids are uncommon. We therefore 



Figure 3. Effective radii of galaxies binned by redshift, for colour 
and Sersic index subsamples. The left column shows low-n galax¬ 
ies, the right column shows high-n galaxies. The top, middle and 
bottom groups of four panels contain red, green and blue samples, 
respectively. Within each group, the top panels show normalised 
histograms, whilst the bottom panels present corresponding cu¬ 
mulative histograms to better visualise subtle shifts in the distri¬ 
butions. In each panel lines in light colours show galaxies in the 
low redshift bin (0 ^ 2 : < 0.15), mid-shades indicate the interme¬ 
diate redshift bin (0.15 ^ z < 0.25), whilst the darkest colours 
show galaxies in the high redshift bin (0.25 ^ z < 0.3). 

consider our measurements of the wavelength dependence 
of effective radius to be robust out to z = 0.3. 

We study in the same manner. In Fig. we show 
the wavelength dependence of Sersic index for galaxies split 
by u — r colour, Ur, and redshift bin. Again, some small 
trends are visible. Galaxies at lower -2 tend to display 
further from unity. This is most significant for low-n galax¬ 
ies, and suggests that V14 may slightly underestimate the 
wavelength dependence of n for these systems. Apart from 
this small effect (approx. 18.5% offset in median values com¬ 
pared to the width of the distributions, as described pre¬ 
viously), we consider our measurements of the wavelength 
dependence of Sersic index to be robust out to 2 = 0.3. 

As in Fig. 17 of V14, we summarise the relationship 
between and TZf for our three redshift bins in Fig. 5 
The small trends with redshift, described above, are again 
evident in this figure. However, for each subsample, the red¬ 
shift bins lie close to one another. In addition to observa¬ 
tional effects, redshift trends may result from real changes in 
the galaxy population. However, given our volume-limited, 
colour-selected samples, and the narrow range of redshift, 
we do not expect evolution of the galaxy population to play 
a significant role in these trends. The distinctions between 
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Figure 4. Sersic indices of galaxies binned by redshift, for colour 
and Sersic index subsamples. Panels and colours are chosen sim¬ 
ilar to FigUl 


subsamples, particularly high- and low-n systems, are main¬ 
tained independent of redshift. The conclusions of V14 are 
thus robust. With that established, we now move on to ex¬ 
amining how the wavelength dependence of galaxy structure 
varies with luminosity. 


3.2 Luminosity-dependence of galaxy structure 

V14 specihcally chose a volume-limited sample of bright 
galaxies with Mr < —21.2. Here we investigate whether the 
conclusions drawn from their study apply to a wider range 
of galaxy luminosities. 

Our sample is divided into three magnitude bins: 
-22.48 ^ Mr < -21.48, -21.48 ^ < -20.48 and 

—20.48 ^ Mr < —19.48. These intervals ensure similar num¬ 
bers of galaxies in each bin, allowing for meaningful compar¬ 
isons. 

In Fig. we present the wavelength dependence of ef¬ 
fective radius for subsamples divided hy u — r colour, Ur and 
luminosity. 

For high-n galaxies, there appears to be a consistent 
trend in the distribution of TZf with changing absolute mag¬ 
nitude. Fainter high-n galaxies have TZf closer to unity than 
brighter high-n galaxies. Their sizes therefore vary less as a 
function of wavelength. The most luminous galaxies, such 
as those studied by V14, display the strongest variation of 
their sizes with wavelength. 

For the low-n samples, although K-S tests indicate some 
differences between the distributions in the luminosity bins, 
there is little obvious consistent variation. The blue sample 


e 

1.0 
^ 0.9 
0.8 
0.7 


0.5 

0.4 

0.2 


Red: > 2.5 


Red: < 2.5 


Green: > 2.5 


Green: n,. < 2.5 


Blue: n,. > 2.5 





-# 




± 


t. 




0.00 < z <0.15 
0.15 < z < 0.25 
0.25 < z < 0.30 


0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Rm/Reia) 


Figure 5. Median versus median for galaxies in each 
of our different subsamples divided by colour, Sersic index and 
redshift bin. Colours are the same as in figure 17 of V14. Error 
bars show uncertainty on the median, estimated as 1.253 (j/'s/]V, 
where a is the standard deviation about the median and N is 
the number of galaxies in the sample. There are small trends 
with redshift due to observational biases. However, the distinction 
between the subsamples, particularly the contrast between low- 
and high-n, remains very clear. Note that the blue 0 ^ z < 0.15, 
n > 2.5 point lies at = 0.15, J\f^ = 0.56, but is consistent 
with the higher z points given its large uncertainties. 


displays similar trends to the high-n galaxies, while typical 
green galaxies show no luminosity dependence. Red galaxies 
hint at an opposing behaviour, such that the most luminous 
bin displays less variation in size with wavelength. 

On average across all samples, the offset of the distri¬ 
butions with respect to the standard deviation of the distri¬ 
butions is 27%, calculated as above, for the redshift distri¬ 
butions. 

Again we study in the same manner; Fig. shows 
the wavelength dependence of Sersic index for the same sub¬ 
samples as Fig.[^ There appears to be a consistent trend in 
the distribution of with changing absolute magnitude 
for both low- and high-n galaxies, with a 47% mean offset of 
the distributions, compared to the width of those distribu¬ 
tions. In the high-n sample, distributions for all luminosity 
bins peak around one, supporting the conclusion of V14 that 
bulge-dominated (and likely one-component) systems show 
little change in Sersic index with wavelength. Comparing 
cumulative distributions for the high-n samples, we see that 
bluer and less luminous samples tend to have slightly lower 
values of . Typical bright high-n galaxies have very 
close to unity. 

The low-n samples consistently display values of 
above one, indicating an increase in Sersic index with wave¬ 
length. This offset from unity strongly increases with lumi¬ 
nosity, such that the brightest low-n galaxies display the 
greatest dependence of Sersic index on wavelength. 

In a similar manner to before. Fig. shows the rela¬ 
tionship between and TZ^ for our three magnitude bins. 
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Figure 6. Effective radii of galaxies binned by absolute magni¬ 
tude, for colour and Sersic index subsamples. Panels and colours 
are chosen similar to FiglU 

High Sersic index galaxies show a clear variation in effective 
radins with luminosity (brighter galaxies showing a greater 
decrease in Rb with wavelength). Low Sersic index galaxies 
show a change in n with wavelength (brighter galaxies show¬ 
ing a greater increase in n with wavelength). The overall ef¬ 
fect is that the differences between low- and high-n galaxies 
become more pronounced with increasing luminosity. 

We have performed all of the analysis in Sec. |3.1| and |3.2| 
using an optical baseline (u-z), although using the same fits 
to the full u-K dataset. As one would expect, over this nar¬ 
rower wavelength range the differences between the samples 
narrow, though not in terms of significance, as the scatter 
shrinks roughly in proportion. All of the behaviour seen in 
the g versus H plots is qualitatively the same, and our in¬ 
ferences would be unchanged. 

In terms of overall colour, concentration and size, early- 
and late-type galaxies become more similar with increasing 
luminosity. However, when one considers the joint distribu¬ 
tion of these three properties, as in this paper, one instead 
finds a divergence which widens with luminosity. The pro¬ 
cess by which late-type (low-n) galaxies grow must promote 
the variation of profile shape with wavelength. The growth of 
early-type (high-n) galaxies must also maintain a radial seg¬ 
regation of their stellar populations, but crucially in terms 
of size, not profile shape. Furthermore, any process which 
transforms late-types to early-types must do so in a man¬ 
ner that erases the wavelength-variation of profile shape, 
but enhances that of size. Realistic mechanisms for galaxy 
growth and transformation must reproduce the trends we 
observe. 


Figure 7. Sersic index of galaxies binned by absolute magnitude, 
for colour and Sersic index subsamples. Panels and colours are 
chosen similar to FiglU 
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R^{H)lR^{g) 

Figure 8. Median J\f^ versus median for galaxies in each 
of our different subsamples divided by colour, Sersic index and 
luminosity. Colours are the same as in figure 17 of V14. 
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8 Kennedy et al. 

3.3 Morphological classifications 


The work of V14 assumes that one can perform a mean¬ 
ingful division of the galaxy population using Sersic index 
as a proxy for structure or morphology. Typically Ur < 
2.5 galaxies are associated with structures with prominent 
discs, whilst Ur > 2.5 galaxies are thought to be spheroid- 
dominated. To verify that our results do not strongly de¬ 
pend on this assumption, we now examine the distributions 
of and TZf for samples selected by visual morphology. 
In Figs. 1^ and we show the wavelength dependence of 
n and for galaxies separated by Hubble types from the 
GAMA Visual Morphology catalogue ( [Kelvin et al.|[2014 |. 
Note that visual morphologies are only available for galaxies 
with 0.025 < 2 < 0.06 and Mr < —17.4. Figurej^shows only 
galaxies within the overlapping redshift and magnitude re¬ 
gions, i.e. 0.025 < z < 0.06 & Mr < —21.2 (V14 sample) or 
Mr < —19.48 (low-a sample). Also plotted (as triangles) are 
the median locations of the high- and low-n populations, 
selected within the same redshift and absolute magnitude 
limits as the morphologically-classified points. Background 
contours show the distribution of galaxies in each of our 
volume-limited samples for which we do not have a visual 
morphological classification. 


Despite the severe limitation in sample size that the vi¬ 
sual morphological classifications impose, they confirm the 
behaviour of our high- and low-n galaxy samples. In both 
volume-limited samples we see that the median location of 
the high-n galaxies lies very close to the median point of the 
elliptical sample, whilst the median of the low-n galaxies lies 
close to that of the Sab-Scd galaxies. Vika et al. (20151 see 
identical behaviour for a sample of very low-redshift galax¬ 
ies, and show that this is robust against artificial redshifting, 
and investigate the use of Af and TZ as morphological clas¬ 
sifiers. 


Note that the offset between the late-type/low-n points 
at low-redshift and the contoured distribution at higher red- 
shifts appears to confirm the (small) biases inferred in Sec¬ 
tion 3.1 At lower-redshift (i.e. in better quality data) we 
measure slightly higher TZ^ and J\f^ for low-n galaxies. Our 
previous results therefore somewhat underestimate the dif¬ 
ference between the two galaxy populations. However, even 
at z = 0.3 the contrast between their behaviour is such that 
our conclusions are unaffected. 


The finer distinctions between the Hubble types reveal 
some further details in the AJ'^-TZf figures. Ellipticals and 
Sab-cd galaxies form the opposite extrema of the distribu¬ 
tions. SO-Sa galaxies are intermediate between the two, but 
closer to the ellipticals. Perhaps surprisingly, late-type disks, 
Sd-Irr, occupy the same intermediate region. This may be 
a result of their similar one-component natures, with higher 
values of AJ"^ indicative of system with comparable bulge 
and disc components. Note that TZ and Af may be noisier 
for systems that depart greatly from an elliptical Sersic, and 
more difficult to interpret physically. Nevertheless, measur¬ 
ing these properties still gives insight into their structure 
and allows us to make qualitative comparisons. So, in Fig. 
[^we see that, on average, Sd-Irr galaxies display proper¬ 
ties similar to late-spirals, but with less peakiness at longer 
wavelengths. It is also worth noting that the presence of a 
bar does not appear to greatly alter the measured wave¬ 
length dependence of galaxy structure. 


10.0 


All galaxies: z <0.3, <-21.2 

Morphologically classified, 0.025 <z <10.06, Mr<-21.2 

Ur classified galaxies, 0.025 <z<0.06, il^r<-21!2 
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Figure 9. AA^ versus TZ^ for galaxies in our high-luminosity 
volume-limited sample with Hubble type classifications in the 
GAMA visual morphology catalogue. Median points for each 
Hubble type bin are plotted as stars, with associated error bars 
(estimated as 1.253iT/\/fV). Morphological classifications are only 
available for a low-redshift subset of galaxies. For comparison, the 
median values for galaxy samples with the same redshift limits, 
but separated by Ur, are plotted as triangles. The distribution of 
the full volume-limited sample (to the same luminosity limit, but 
extending out further in redshift) is plotted as contours. 



Figure 10. AA^ versus TZ^ for galaxies in our lower-luminosity 
volume-limited sample. Points and shading are as described for 


Fig.[9| 
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3.4 The effect of dust 


Late-type galaxies contain significant quantities of dust, 
which can strongly affect their measured properties (e.g., 


Popescu et al.|2002||Pierini et al.|2004||Holwerda et al.|2005l 


2009 Mollenhoff, Popescu & Tuffs 20061. The effects vary 


with the amount and distribution of dust within the galaxy, 
as well as the observed inclination and wavelength. 

A wide range of observations, together with detailed 
modelling ( [Popescu et al.pOOOp , have helped to establish a 
typical geometry for the distribution of dust in the disks of 
star-forming galaxies ( Tuffs et al.||2004 l. When accounting 
for attenuation, this dust model can be parametrised by a 
single dominant parameter: the central face-on optical depth 
in the B-band, r, while remaining widely applicable, at least 
to relatively massive spiral galaxies. It should, however, be 
noted that in lower-mass galaxies the dust distribution may 
be significantly different ( Holwerda et al.||2012 Hinz et al. 
|2007| ). The [Tuffs et al.| ( |2004p model has been used to quan¬ 
tify, and potentially correct for, the effect of dust on galaxy 


structural parameters (Mollenhoff, Popescu & Tuffs 2006 
[Pastrav et al.|2013a|b[). 


Figure 


11 


shows the effect of dust on M? and TZtf for 


pure discs fit with a single-Sersic profile, as determined by 


Pastrav et al. (201 Sal. Their results are shown for a range 


of inclinations along loci corresponding to a wide variety of 
central dust opacities. Typical spiral galaxies have r ~ 2-4 


( Keel fc White Ill|[2001 Holwerda et al.[[2005 Driver et al. 


2007b Masters et al. 20101. Projection effects, as defined 
and predicted in Pastrav et al. ( 2013a|b I, are not considered 
here, as their contribution is minimal compared to dust ef¬ 
fects. As the purpose of this figure is illustrative, we have 
simply taken the published K and B-band dust effects from 


Pastrav et al. (2013al, without attempting the small inter¬ 


polation to the g- and B-bands. These dust effects are shown 
as applied to a nominal ‘galaxy’ with {JV^, Tiff) = (1, 1), 
and for each r value we show how & TZf change with 
inclination. Fig. shows that dust always acts to decrease 
TZf, i.e. increase apparent size more at shorter wavelengths. 
It also tends to increase , i.e. raise the Sersic index with 
increasing wavelength, except for at high opacities and close 
to edge-on inclinations. 


Note that the stellar disk emitting in the optical bands 
in the [Tuffs et al.[ ( [2004[ ) model actually has an intrinsic 
stellar population gradient, such that the scale length of 
the old stellar disc decreases by ~30 percent from the B 
to K bands, corresponding to TZf ~ 0.7. Including this in¬ 
trinsic TZ^ would make the loci overlap with the observed 
distribution in Fig. [m One can see that for a population of 
pure discs with a variety of opacities and inclinations, a cor¬ 
relation in would arise from their variation with 

inclination, while a scatter about that correlation would 
be associated with variations in opacity. However, as one 
would expect, dust effects alone cannot account for the full 
observed distribution. Late-types extend to sub¬ 

stantially greater than can be attributed to dust, pre¬ 
sumably as a result of the presence of a central bulge that 
becomes more dominant at longer wavelengths. However, 
there remains the possibility that more varied dust models 
may have more success in reproducing our observed depen¬ 
dence of A/* and IZ. For example, the dust within lower-mass 
galaxies is more likely to be vertically distributed through- 
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Figure 11. Dust corrected J\fR vs for discs. Underlayed are 
contours for all galaxies in the V14 sample, as in Fig.[^ Overplot¬ 
ted are the mean locations of low- and high-n galaxies within the 
V14 volume limited sample. Black lines show model tracks for op¬ 
tically thin discs having a nominal value (A/*^,7^^)=(l,l), whilst 
coloured lines correspond to model tracks for optically thick discs 
The intensity of the colour in each line shows the inclination of 
the galaxy. Crosses indicate the point at which (1 — cosi) = 0.5. 





out the stellar disk, rather than concentrated in the central 
plane (Holwerda et al. 20121. On the other hand, early-types, 
which generally contain little dust, must have stronger in¬ 
trinsic stellar population gradients than those in the discs. 


4 DISCUSSION 

In this paper we have studied how the wavelength depen¬ 
dence (from restframe u- to H- bands) of galaxy struc¬ 
ture varies for samples selected in a variety of ways. As 
in V14, galaxy properties have been measured using Meg- 
aMorph techniques to fit consistent, wavelength-dependent, 
two-dimensional Sersic profiles in multiple wavebands simul¬ 
taneously. Our results are summarized in terms of the frac¬ 
tional variation in Sersic index and effective radius between 
the g- and B-bands, which we denote and Tiff ■ 


4.1 The physical meaning of TZ and N 

Before discussing the implications of our results, and com¬ 
paring to other studies, we hrst review the meaning of these 
quantities in terms of the appearance and physical structure 
of a galaxy. A value around unity, for either quantity, indi¬ 
cates that that particular aspect of structure does not vary 
with wavelength. TZ gives the variation in galaxy size (Be) 
with wavelength. The vast majority of galaxies are found 
to have TZ < 1, indicating that they are larger in the blue. 
These galaxies therefore display redder colours in their cen¬ 
tres, and become bluer at larger radii. In the rarer case, TZ 
> 1, galaxies are larger in the red, and hence bluer in the 
centre, becoming redder with increasing radius. 

M indicates how a galaxy’s proHle shape (n), or equiv¬ 
alently its central concentration, depends on wavelength. 
Galaxies with Af > 1 have a high-n proHle in redder bands. 
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10 Kennedy et al. 


They are therefore ‘peakier’ in their centres and have more 
signihcant surface brightness at large radii. Galaxies with 
A/” > 1 are thus expected to be redder in the centre and in 
their outskirts, and bluer at intermediate radii (~ Ra). In 
practice, the outermost red regions may be too faint to been 
seen in normal depth imaging, and so the visible outskirts of 
such a galaxy may appear blue. However, the difference with 
respect to 7?. < 1 is that the colour gradient should become 
shallower with radius. Galaxies with Af < 1 obviously show 
the opposite trend, appearing more strongly peaked at bluer 
wavelengths. 

Of course, both Af and TZ may both depart from unity, 
indicating variations in both size and profile shape with 
wavelength. In that case, the resulting colour profile is a 
combination of the two behaviours. These may, to a degree, 
counteract (if Af and TZ are correlated) or reinforce (if anti¬ 
correlated) one another. Indeed, some of the correlation in 
the Af-TZ plane, within subsets of the galaxy population, 
may be a result of this degeneracy. 

We study the effect of re dshi ft on our ability to mea¬ 
sure Af? and TZ^ in Section 


NIR-optical colour gradients. This agrees with our Ending 
that Rc depends on wavelength more strongly for more lu¬ 


3.1 


For low-n (i.e. typically 
late-type) galaxies we hnd a bias of TZf and Af? to lower 
values {TZf away from unity, Af? toward unity) at higher 
redshifts. However, these observed biases are significantly 
smaller than the differences between our galaxy subsamples, 
particularly for low- versus high-n galaxies. The conclusions 
of V14 therefore are not significantly affected by the gener¬ 
ous redshift limit adopted in that study. The quantitative 
results of V14 are also relatively robust and may be com¬ 
pared to other studies of comparable galaxy populations in 
similar quality imaging. However, when comparing to much 
lower redshift samples of late-types, e.g. in Fig. |10[ the red¬ 
shift bias must be borne in mind. 

In Section |3.2[ the wavelength-dependence of both R^ 
and n are found to depend on luminosity, with brighter 
galaxies generally showing stronger trends with increasing 
wavelength: to smaller sizes (especially for early-type galax¬ 
ies) and higher-n (especially late-type galaxies). The con¬ 
trast between the behaviour of high- and low-n galaxies 
lessens with decreasing luminosity, suggesting greater struc¬ 
tural similarity between (Sersic index selected) early- and 
late-type galaxies at lower luminosities. 


4.2 Comparison of observed trends of R and N 
with other studies 

Our work is consistent with that of ILa Barbara et al.l 


12010b I. We find an extremely similar decrease in effective 


radius for red, n > 2.5 galaxies at the same luminosities. 
Furthermore, we have shown that this is robust to redshift 
effects and the use of visual morphology, rather than Sersic 
index, to select early-types. Our results also agree with those 
of |Kelvin et al.| ( |2012[ ) and |Ko fc Im| ( |2005[ ), both of which 
find a nearly 40 per cent decrease in Rc over the correspond¬ 
ing wavelength range. 

The existence of negative colour gradients in elliptical 
galaxies, particularly strong in optical-NIR colours, has been 


known for a long time (e.g., Peletier, Valentijn & Jameson 
p^ .[La Barbera et al.| ( |2010a[ ) use their Sersic models, in- 
dependently fit to each wavelength, to determine colour gra¬ 
dients for their sample of early-type galaxies. They find that 
more optically-luminous galaxies display stronger (negative) 


minous high-n galaxies. However, La Barbera et al. (2010a I 


see little variation in colour gradients with NIR luminosity 
or stellar mass. This is consistent with the bluer stellar pop¬ 
ulation being located at larger radii, and this feature being 
more prominent in more luminous galaxies. 


La Barbera et al. (20121 show that, compared to optical- 


NIR measurements, colour gradients based only on optical 
bands are weaker. They are nevertheless widely observed 


(e.g., Gonzalez-Perez, Gastander & Kauffmann 20111. Op¬ 
tical gradients also do not show as much variation with lu¬ 
minosity, being strongest for intermediate luminosities (La 
Barbera et al. 2012| Roche, Bernard! & Hyde 20101. Us¬ 
ing stacked optical (SDSS) images, D’Souza et al. (20141 
find strong evidence for the presence of blue (in g — r) 
halo components in various galaxy populations. This ex¬ 
tended component can account for around half of the light 
of a galaxy, being more prominent in higher-n and more lu¬ 
minous galaxies. Even the majority of massive early-type 
galaxies at z ~ 1.5 are found to display negative colour gra¬ 
dients ( Gargiulo et al.|[2012 (. 

A small subset of local early-type galaxies display pos¬ 
itive optical colour gradients, due to the presence of blue 


cores associated with recent star-formation (Suh et al. 20101. 


These become more prevalent at higher redshifts (Ferreras 


et al. 20051. Such galaxies correspond to our high-n, very 


blue, selection. These galaxies display TZg < 1, like other 
early-types, implying they are more compact in the red. 
However, in contrast to our other subsamples, this is com¬ 
bined with Af? < 1, implying their profiles are more 
strongly peaked in the blue. The intermediate and outer gra¬ 
dients of these galaxies are therefore apparently like other 
ellipticals, but they must display a blue excess in their cores. 

Typical spirals (Sab-Scd) show less Rc variation with 
wavelength, compared to early-types, but it is still signif¬ 
icant: a ~ 20 per cent decrease from g to H. This value 
has little dependence on galaxy luminosity, at least over the 
range we probe here. However, note that the wavelength 
variation of Ra for late-types is somewhat susceptible to 
being overestimated in lower-quality data. Both early-discs 
(SO-Sa) and very late-types (Sd-Irr) display behaviour inter¬ 
mediate between ellipticals and spirals. 

In terms of Sersic index, we find that early-type (high-n) 
galaxies show little variation in n, while late-types substan¬ 
tially increase in n with increasing wavelength. Again, early- 
discs (SO-Sa) and very late-types (Sd-Irr) fall between the 
two behaviours. Our results for disc galaxies agree well with 


those of Kelvin et al. (20121, who find that Sersic indices 


almost double from u to K, although most of the change 
occurs over the optical range. 

Various studies have found the Sersic indices of early- 
type galaxies to be generally unchanged with wavelength 
(Kelvin et al. 2012 La Barbera et al. 2010a| 20121, per¬ 
haps with a slight increase in n with wavelength. The esti¬ 
mated changes in n and Ra with wavelength for these, and 
other, studies are displayed in Table We do not generally 
see such an increase, rather a small decrease is more com¬ 
mon in our early-type samples. However, redshift appears 
to have an effect, as does luminosity. In our lowest-redshift, 
brightest samples we find Sersic indices that are constant, or 
slightly rising, with wavelength. Our visually-classified ellip- 
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‘Early type’ 

‘Late type’ 

Notes 


n 

Re 

n 

Re 


LBIO 

i*t 

28 4 

- 

- 

5080 bright spheroids from SPIDER, g- to X-band. ||La Barbera et aL|2010a|l 

Kell2 

30 t 

38 4 

52 t 

25 4 

167,600 galaxies from GAMA, g- to R-band. l Kelvin et al.|2012l 

V14 

l*t 

45 4 

38’’'t 

25 4 

14,274 galaxies from GAMA-II GOO region, Mr < —21.2,2: < 0.3, u- to i7-band. 

L15: 10® Mq 

- 

- 

- 

16 4 

8399 galaxies from GAMA-II, 0.01 < 2 : < 0.1, g- to Rs-band. 1 Lange et al.|2015l 

L15: IOIOMq 

- 

13 4 

- 

13 4 

As above. 

L15: 1O“AL0 

- 

11 4 

- 

- 

As above. 

K15: Mr ~ -20 

12 4 

23 4 

29 t 

13 4 

5331 galaxies from GAMA-II GOO region, 2 < 0.15, g- to i7-band (this study). 

K15: Mr ~ -21 

5t 

25 4 

40 t 

15 4 

As above 

K15: Mr ~ -22 

5t 

33 4 

55 t 

12 4 

As above 


Table 3. Comparison of results from previous studies. Each number shows the percentage change in n and iJe between the wavebands 
given under ‘notes’, f = increase; = decrease; ““ = value has been calculated from information given in the paper; - = no data available... 


ticals also have very close to unity. On the other hand, 
fainter and bluer high-n galaxies tend to display Sersic in¬ 
dices that decrease with wavelength, by several tens of per 
cent. La Barbera et ah] (20121 find that larger (and hence 


typically more luminous) early-type galaxies show stronger 
trends of n increasing wavelength, matching the form of be¬ 
haviour we see. 

The redshift biases and luminosity trends we have 
found allow us to, at least partly, reconcile the dramatic 
wavelength-dependence reported by V14 with the more 
modest behaviour seen by some other recent studies. |Kelvin| 


et al. (20121 find slightly weaker trends than V14, but they 


use a magnitude-limited sample, dominated by somewhat 
less luminous galaxies. When we adopt the same selection, 
our results are in excellent agreement. 


Lange et al. (20151 find a substantially weaker depen¬ 


dence of on wavelength, though it is still signihcant. We 
have shown that dividing galaxies by visual morphology (as 


Kelvin et al. |2014|) confirms the behaviour inferred from 


selections based on colour and Sersic index, but severely lim¬ 
its the useable sample size. 

For late-type galaxies, our results are actually in good 
agreement with Lange et al. (20151. If we consider our low 
redshift sample, we find a similar ~ 15 per cent decrease in 
size. At higher redshifts we appear to slightly overestimate 
the strength of the size decrease. However, for early-types, 
there is a larger discrepancy, although the qualitative be¬ 
haviour is the same. Redshift (i.e. data-quality) appears to 
have little effect on our results. Considering our Mr ~ —20 
(log M* ~ 10) selection we still find a ~ 25 per cent decrease 
in the effective radii of red high-n galaxies. This is not re¬ 
duced if we limit our analysis to visually-classified ellipticals. 
The wavelength-dependence results featured by|Lange et al.| 
( [Ms] ) are based on linear fits to the Rc-Mr, relation, and 
hence galaxies with relatively low masses, logM* ~ 10, ex¬ 
ert a strong influence. These different analyses may have 
something to do with our disagreement. However, it may 
also be that the difference in depth between the UKIDSS 
LAS and VIKING VISTA NIR imaging is responsible, as 
suggested by [Lange et al.| (|2015|). Many, but not all, of the 


studies that also find strong trends for early-type galaxies, 
as mentioned above, are also based on SDSS and UKIDSS 
LAS data and hence may be similarly affected. In the near 
future we will have the results of our methodology applied 
to the same SDSS plus VIKING data set, so will be better 
placed to evaluate this issue. 


4.3 Connection to stellar population gradients 


The wavelength dependence of structure that we have mea¬ 
sured, i.e. A/” and TZ, must be driven by radial variations 
in stellar populations and/or dust extinction. For early-type 
galaxies, many studies have found that colour gradients are 
caused by negative metallicity gradients. Age is generally 
fairly constant with radius, or even slightly increasing, and 
hence acting against the metallicity trend. The outskirts of 
massive early-type galaxies (ETGs) are therefore typically 
more metal poor and older than their cores, with weaker gra¬ 


dients for lower mass ETGs. (e.g.. La Barbera et al. 2012 
and references therein). 

Metallicity gradients are expected in models of mono¬ 


lithic collapse (Worthey, Trager & Faber 19951, and this 


mechanism may be responsible for the gradients in the cen¬ 
tres of today’s galaxies. However, most massive galaxies are 
expected to have experienced a major merger since their for¬ 
mation, which will have partly erased these initial gradients. 
It has been discovered that early-type galaxies must at least 
double in effective radius between z ^ 2 and today (e.g., 
see [van der Wei et al.||2014| . Although various processes 
have been proposed, the dominant mechanism appears to 


be minor mergers (Hopkins et al. 2009 Lackner & Gunn 


20121. The stellar material accreted in such mergers tends 


to contribute most significantly to the outskirts of the more 


massive galaxy, building up an outer envelope (Huang et al. 
20131. The effect is to increase effective radius over time and 


produce stellar population gradients. 

Our measured increase in effective radius toward shorter 
wavelengths, which is stronger for more luminous systems, 
{TZf ~ 0.55-0.75 from Mr ~ —22—20), is consistent with 
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observations of negative colour gradients and extended blue 
components Barbera et al.|[2012[ ). The lack of varia¬ 
tion in Sersic index with wavelength (A/” ~ 1) suggests that 
the extended profile of the stellar populations seen in the 
blue is structurally consistent to the profile traced by those 
stellar populations dominating at red wavelengths, possibly 
indicative of a similar origin. Together, these observations 
are consistent with a gradual build-up of massive early-type 
galaxies through accretion of low-mass galaxies. 

Lower-luminosity, and bluer, early-type (high-n) galax¬ 
ies tend to show slightly peakier profiles in the blue, which 
may be an indication of an additional blue stellar compo¬ 
nent in their centres, possibly as a result of recent central 
star-formation. 

Dust in early-type galaxies is generally not expected to 
be present in sufficient quantities to have a significant effect 
on their observed structure. However, Rowlands et al. (20121 
do hnd that ~ 5 per cent of luminous early-type galaxies are 
detected in the infrared, implying that they must contain 
signihcant fractions of dust (also see Bendo et al. 20061. 
It is uncertain how this dust impacts upon their measured 
structural parameters, but it may be responsible for some of 
the scatter we see in A/” and TZ for our early-type samples. 

In late-type galaxies, dust is likely to play a more dom¬ 
inant role. The apparent differences in optical structure in 
disc galaxies across the Hubble sequence is largely deter¬ 
mined by the presence of young stars and dust. The old 
stellar discs of spiral galaxies, as seen in the NIR, are very 
similar at all Hubble types ( Block et al.||1999 |. 

Modern studies of the opacity of spiral discs find 
that they suffer from substantial face-on and inclination- 
dependent extinction (e.g., Holwerda et al.| 2005 Masters 
et aI.|[20T0|. Testing the model of 


Tuffs et al. (20041 with a 


large galaxy sample. Driver et al. (2007b I find good agree¬ 
ment and infer that bulges and the central regions of discs 
suffer from substantial attenuation of their optical light. 
This has a signihcant impact on measurements of their struc¬ 


tural parameters (Mollenhoff, Popescu & Tuffs 2006 Pastrav 
|et al.|[2013a|b| ). However, as we have shown in Section |3.4| 
dust alone cannot explain the observed wavelength trends. 
It seems obvious that the two-component, bulge-disk, nature 
of late-type galaxies, with their contrasting stellar popula¬ 
tions and structures, must have at least some impact on the 
wavelength-dependence of their overall structure. We will 
investigate this further in a forthcoming paper. 


5 SUMMARY 

The variation of galaxy structure with wavelength reveals 
the connections between the stellar populations within a 
galaxy and their spatial distributions. Specihcally, we con¬ 
sider the change in Sersic index and effective radius from 
g- to iL-band, which we denote A/” and TZ. As shown by 
V14, the majority of early-type (high-n) galaxies show lit¬ 
tle variation in their Sersic index with wavelength {Af ~ 1), 
but are significantly smaller at longer wavelengths {TZ < 1). 
This behaviour is suggestive of a structure formed by simi¬ 
lar (violent) mechanisms being apparent at all wavelengths, 
but with a scale that strongly depends on the colour of the 
stellar population. On the other hand, late-type galaxies 
(low-n) display a substantial increase in Sersic index with 


wavelength (A/” > 1), suggesting a variation in the type of 
structure dominating the light at different wavelengths, but 
less variation in size {TZ < 1). Very blue, high-n, galaxies 
present an interesting contrast, having signihcantly higher 
Sersic indices in the blue, perhaps indicating the presence of 
a peaky blue component. 

In this paper we have followed V14, using optical- 


NIR imaging and redshifts from the GAMA survey (Driver 
et al.|2009 (, with multi-band single-Sersic fits performed us¬ 


ing tools developed by the MegaMorph project (Bamford 


et al. 20121, in order to study four further aspects of the 


wavelength-dependence of galaxy structure. We consider the 
distributions of N and TZ for a variety of volume-limited 
samples subdivided by u — r colour, Sersic index and mor¬ 
phology. The main conclusions from our analysis are: 


• We have verihed that our measurements of N and TZ 
are robust to the effects of redshift, strengthening our earlier 
results (V14). Early- and and late-type galaxies (selected us¬ 
ing Sersic index) present strongly contrasting behaviour in 
terms of the wavelength-dependence of their structure. Out 
to 2 : ~ 0.3 we see no substantial changes in TZ or A/”, par¬ 
ticularly when compared to the striking differences between 
galaxy subsamples. What small biases are present, suggest 
that at higher redshift (lower S/N and poorer resolution) we 
may be slightly overestimating the offset of TZ from unity, 
and underestimating the offset of N from unity, for low-n 
galaxies. Both of these effects act to reduce the apparent 
difference between the high- and low-n populations. The es¬ 
timated contrast is therefore even more pronounced than 
that determined by V14. 

• The strengths of A/” and TZ depend on galaxy luminosity: 


- High-n galaxies with lower luminosities have TZ closer 
to unity than their brighter counterparts. Most high-n 
galaxies show a peak in M around unity for all lumi¬ 
nosities, supporting the suggestion of V14 that these are 
single-component objects. 

- Low-n galaxies display weaker, mixed trends in TZ 
than high-n galaxies. There is a striking trend in the N 
distributions of low-n galaxies with luminosity: brighter 
objects display M further from unity, for all colour sub¬ 
samples. 


• The interpretations of V14 are supported by visual 
morphological classifications of a low-redshift subset of our 
sample from Kelvin et al. (20141; our low-n samples share 
the same part of the Af~TZ diagram as Sab/Scd galaxies, 
whilst the high-n samples follow the distribution of ellipti¬ 
cal galaxies. Both early- and late- disc galaxies occupy in¬ 
termediate regions of the N-TZ plane, confirming that the 
extreme values for intermediate types are related to their 
two-component nature. 


• Results from fitting dusty galaxy models (Pastrav et al. 


2013a I indicate that some of the wavelength-dependence of 


disc galaxy structure may be attributable to dust. The nat¬ 
ural distribution of disc inclinations may account for the 
trend seen in the N-TZ plane, while varying central dust 
opacities may account for the scatter in this relation. How¬ 
ever, pure dusty discs cannot reach the values of N observed 
for intermediate-type spirals. 


Further improving our understanding of the 
wavelength-dependence of galaxy structure will require 
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more detailed consideration of the two-component nature 
of galaxies. We are currently exploring this topic using 
bulge-disc decompositions, and these results will be the 
subject of Kennedy et al. (in prep.). 
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